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Abstract —In this letter, we proposed a shared-antenna full- 
duplex framework for the multiuser MIMO system with a large- 
scale antenna array exploited at the base station to transmit 
and receive signals simultaneously. It has the merits of both 
the full-duplex system and the time-division duplex massive 
MIMO system, i.e., the high spectral efficiency and the channel 
reciprocity. We focus on the zero-forcing (ZF) and the maximal- 
ratio transmission/maximal-ratio combining (MRT/MRC) linear 
processing methods, which are commonly used in the massive 
MIMO system. As the main finding, we prove that the self¬ 
interference in the completely dependent uplink and downlink 
channels of this full-duplex system can be suppressed by the 
large-scale antenna array. 

Index Terms —Shared-antenna, in-band full-duplex, MU- 
MIMO, self-interference, massive MIMO. 

I. Introduction 

ARGE-SCALE MIMO (massive MIMO) and in-band 
full-duplex have been attracting lots of interests from 
academia and industry nowdays m-tm because both of 
them can achieve high spectral efficiency. In a traditional full- 
duplex system 0-0, a separate antenna or antenna array is 
employed to receive the signals from the remote transmitters, 
meantime the self-interference (SI) is also received from its 
own transmitter nearby. Compared with the desired signals, 
the SI is much stronger, thereby it must be suppressed before 
decoding the signals from the remote transmitters. Many 
studies have shown that the SI can be mitigated by active 
or passive suppressions, in the wireless propagation, analog 
circuit, or digital domain 0. 0-0. Meanwhile, a massive 
MIMO system usually works in the half-duplex mode, such as 
the time-division duplex (TDD) and employs a single antenna 
array so that the channel reciprocity can be assumed 0 , a, 
0. In a TDD massive MIMO system, the downlink channel 
matrix is a transpose of the uplink channel matrix, therefore, 
the channel reciprocity can help reduce the burden of channel 
state information (CSI) feedback. With a large-scale antenna 
processing method, the interference and noise will be vanished 
when the scale of antenna array is large enough 0. The 
combination of the full-duplex and the massive MIMO is 
actually becoming a hot topic recently j9j. 
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However, this combination faces the conflicts between the 
separate antenna arrays in a full-duplex system and the single 
antenna array in a massive MIMO system. Two types of full- 
duplex antenna interfacing methods have been summarized 
in 0, namely, the separate-antenna and the shared-antenna 
interfacing. The shared-antenna interfacing can reduce the 
antenna pair to a single antenna with a circulator, which 
performs the directional routing for the electromagnetic wave 
in a full-duplex wireless transceiver. The circulator has been 
also introduced in flOl . In this letter, we propose a framework 
for the shared-antenna full-duplex massive multiuser MIMO 
system, it has a single antenna array equipped with circulators. 
This model has the merits of both the full-duplex and the 
TDD massive MIMO, i.e., the high spectral efficiency and 
the channel reciprocity. When the antenna array scale M is 
large enough, we find that the SI in a shared-antenna full- 
duplex massive MIMO system can be suppressed by the linear 
processing methods regardless of the complete dependence of 
the uplink and the downlink channels. The simulation results 
are provided to verify the effects of our proposed model in 
suppressing the SI. 

Notation: The superscripts [•]*, [-] T , and [-] ff stand for the 
conjugate, transpose, and Hermitian transpose of matrices or 
vectors, respectively. The fr(-) represents the trace of a square 
matrix. The expectation operator is denoted as E{ ■}. The [-}ij 
represents the (*, j)th entry of a matrix. The represents 
almost sure convergence, and u n = 0(w n ) denotes that there 
exists a constant P, such that u n < Pw n for all n. 

II. System Model 

A. Shared-Antenna Full-Duplex Massive MU-MIMO Model 

Fig. 1 shows the shared-antenna full-duplex massive MU- 
MIMO system, where K user equipments (UEs) are served by 
a base station (BS). Each UE is equipped with one antenna, 
while the BS has M antennas. Both the UEs and the BS 
are working in the full-duplex mode, i.e., they can transmit 
and receive signals using the same time-frequency resource. 
Like the separate-antenna full-duplex system 0, the physical 
isolation is exploited between the BS antennas to avoid the 
front end saturation at the receiver. Therefore, the Mx 1 uplink 
received signal vector at the BS is given as 

Vbs = VKGx u +^p^GsS + n, ( 1 ) 

where G is the M x K channel matrix between the K 
users and the BS; x u is the vector of symbols transmitted 
by the K users, and it has unit variance E{at u a:^} = Ik', 
The normalized uplink power is p u ; G s is the M x M SI 
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Fig. 1. The shared-antenna full-duplex massive MU-MIMO model. 


channel matrix between the BS transceivers; s is the downlink 
precoded vector transmitted to the K users, and it has a 
total expected power E{s fl s} = 1; The normalized downlink 
power is pd', n is an additive complex Gaussian noise vector 
with zero mean and unit variance. 

The channel matrix G incorporates the small-scale Rayleigh 
fading and the large-scale fading, so we have 

G = HD 1 ? 2 , (2) 

where H is the M x K matrix with i.i.d. entries, which 
represent the small-scale fading between the I\ users and the 
BS, and the entry h m k is a circularly symmetric complex 
Gaussian (CSCG) random variable with zero mean and unit 
variance, i.e., h m k ~ CA/"(0,1); D is a I\ x K diagonal matrix 
of the large-scale fading, where [D]kk = Pk- 

The SI channel matrix G s models the direct-path interfer¬ 
ence and the reflected-path interference factors. Unlike the 
single shared-antenna, for the shared-antenna array, the direct- 
path SI comes from not only the circulator leakage, but also 
the antennas inter-coupling if they are not ideally isolated; The 
reflected-path SI is from not only its own reflection, but also 
the others. We assume that the direct-path SI is deterministic 
and the reflected-path SI is random. Therefore, the SI channel 
matrix can be modeled as 

G s = G s + G s , (3) 

where G s is an M x M complex deterministic matrix with 
[Gsjij = Cij, which represents the antennas intra-leakage 
and inter-coupling coefficient (direct-path). G s is an M x M 
complex random matrix (reflected-path) ED, 10. Like the 
signal channel matrix G, G s can be expressed as 

G s = H s Dl /2 , (4) 

where H s is an M x M matrix with i.i.d. entries, which rep¬ 
resents the SI small-scale fading between the BS transceivers, 
it is caused by the reflection from the scatterers. We suppose 
[H s ]ij = hij, where h, :/ is a CSCG random variable with zero 
mean and unit variance; D s is an M x M diagonal matrix 
of the large-scale fading coefficients, where [D s ]mm = 


We assume that all the SI large-scale fading coefficients have 
the same value /?, because the size of antenna array is much 
smaller than the distance of the scatterers from the array. 


B. Uplink Decoding Under the Downlink Precoding SI 


As mentioned above, the downlink precoding can be used in 
the data transmission phase so that the strength of radio signal 
at the corresponding UE can be enhanced. However, the CSI 
should be available at the transmitter before the precoding. 
In this letter, we suppose the BS knows the CSI perfectly 
through a training process. The BS transmitters are assumed to 
adopt the ZF or the MRT precoding scheme for the downlink 
transmission. Correspondingly, the BS receivers use the ZF 
or the MRC linear processing method for the uplink data 
reception. 

Suppose Xd is the K x 1 data vector sent to the users by 
the BS, and s = Ax d, where A is the precoding matrix 0, 

f a ZF G*(G T G*)-\ for ZF, 

A = \ \ (5) 

I £*mrtG , for MRT, 


and Ejajd^d } = The normalization factors azv and 

ctMRT satisfy E{s ff s} = 1, so we have azF — J , 

V = 1 “k 

and «mrt — yj tr?r gh Then, Eq. Q} can be 


rewritten as 


Vbs = VKGx u + ^/p^GsAxd + n. (6) 


In order to obtain the uplink data sent from the users, the BS 
receiver will multiply a linear matrix W T before the received 
signal, which is a function of the signal channel 0. The linear 
receiver matrix can be given as 


W T 


(G h G)~ 1 G h , for ZF, 

G". for MRC. 


(7) 


With the linear receiver at the BS, the uplink signal trans¬ 
mitted from the UEs can be decoded from the received signal 
j/bs- The uplink decoded signal can be expressed as 


r = W T y B s 

= Vp^W t Gx u + ^W t G s Ax d + W T n. (8) 


Correspondingly, the uplink decoded signal for the fcth user, 
namely, the fcth element of r, can be expressed as 

K 

rk = y^KwlgkX^k + X] w k9iXu,i 

i=l,i^k 

+ y/pdwl G s Ax d +wln, (9) 

where and gi, are the fcth columns of W 1 and G, 
respectively; a; u ,fc is the fcth element of a: u ; ^/p^w^gkX^k 
is the desired signal from the fcth user; y/p^^^k 
is the interference from other users; y/pdW^G s Axd is the SI 
and w^n is the additive white noise. 


C. Downlink Decoding 

Like the uplink, there exists SI in the downlink of the full- 
duplex system. The downlink SI includes the inter-user SI and 
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the intra-user SI. The inter-user SI occurs when the user is 
receiving the signal from the BS while other users nearby 
are transmitting signals, so it comes from the reflection and 
the direct coupling of the others’ transmit signals. One option 
for the avoidance of the strong inter-user SI is to switch the 
duplex mode if the users are located in the same place. The 
corresponding UEs can switch their full-duplex mode to the 
half-duplex mode, which are scheduled by the BS. In this 
letter, we assume that all the users are located in different 
places and all the UEs are working in the full-duplex mode, 
and their inter-user SI are attenuated by the distance so that 
no front-end saturation happens. The intra-user SI is from the 
circulator leakage and its own signal reflection, and we assume 
that all the UEs have the intra-SI as well. Therefore, the K x 1 
downlink received signal vector is given as 

VUE = y/PdG T Ax d + y/p^G' s x u + n d , (10) 

where G' is a K x K matrix, and the entry [G'] pg is 
(c' pq + (3 pq h pq ), where c' pq is the intra-leakage and inter¬ 
coupling coefficient (direct-path), (3 pq is the large-scale fading 
of the reflected-path SI, h pq is the small-scale fading, and 
h pq ~ CJ\f{ 0,1); rid is the downlink additive complex 
Gaussian noise vector with zero mean and unit variance. 

By substituting Eq. © into Eq. ( ITOb . and by using the 
asymptotic orthogonality property in EQ, namely, G Ef —‘—A 
D, we yield the kth user received signal with ZF and MRC 
methods 

K 

Vk = \JPdPk%d,k + ^2( C 'kq + fikqh'kq) X u,q + R-d./o (H) 

9=1 

where pk is the linear processing factor, it is given as 


Pk 


ctzf, for ZF, 

auKrMPk, for MRT/MRC. 


( 12 ) 


Note that the desired signal <Jp7\f)kX<\.k in Eq. (fTIT) grows 
with the increase of M while the variance of the downlink SI 
and the noise is bounded. 


III. Self-Interference Suppression 
A. SI Suppression with a Large-Scale Shared-Antenna Array 

In the traditional large-scale antenna array model, the 
interferences are asymptotically orthogonal to the subspace 
spanned by the desired signals when M approaches infinity. 
Therefore, the desired signal can be recovered by projecting 
the received signal in a desired signal subspace with linear 
processing methods such as ZF or MRT/MRC (8j, [0. 

However, in a shared-antenna full-duplex massive MU- 
MIMO system, the interferences are different from those in 
the half-duplex or the separate-antenna full-duplex MIMO 
system. This is because the receivers are very close to the 
transmitters and the antenna array scale is large, the SI is 
much stronger than the other interferences such as the inter¬ 
user interference in the TDD half-duplex system. Moreover, 
the uplink channel is completely dependent with the downlink 
channel because the transmitter and receiver are using the same 
antenna array, which is different from the separate-antenna 
full-duplex system such as 0- Even so, the further study 


shows that the asymptotic orthogonality still holds true in this 
new model, namely, the SI can be suppressed by the projection. 

Lemma 1: Let B denote an M x M complex matrix which 
satisfies tr(BB H ) = 0(M 2 ); Let x = [xi, ..., xm] t and 
y = [y l, Vm] T be M x 1 complex Gaussian random vectors 
of i.i.d entries which have zero means and unit variances. The 
entries of x , y are mutually independent, and they are all 
independent of B. With M approaching infinity, we have 

WT2 xHBx * ^ 0, (13) 

-^x H By ^ 0. (14) 

Proof: Assume x [-> and x,\ are the real part and imaginary 
part of x. They are mutually independent and have the same 
distribution. Consider A/ j /2 x n Bx* , we get 


1 

M 3 / 2 


x h Bx* 


Jf3j5(x£Bx K -x?Bx 1 
- jx^Bx i - jxi Bxr). 


(15) 


Like the proof of Lemma 13 in sm, we have 


E 


x^Bx R - 


M 3/2 


trB 


C 4 


< j^(trBB"y(vj + v 8 ) = 0(1/M 2 ), 


(16) 


where C,\ is a constant, v,\ and v 8 are the 4th and 8th moments 
for the entries of xr. Therefore, the right-hand-sdie (RHS) 
of (fl6l > is summable. By Lemma 3 in ||121 . it follows that 
a/ 3/2 x r Bxr — M l/ 2 trB —A 0. The same result can be 
obtained for xj. Like the proof of Theorem 3.7 in fl3l , we 
have M \/ 2 Xr Bx\ 0. Substitute all the results into Eq. 
we have ( fl3l . Similarly, we can prove ( 1 1 4[ i. ■ 

Remark: The results hold true for either the deterministic 
or the random matrix B. 

Theorem 1: In the shared-antenna full-duplex massive MU- 
MIMO system, the direct-path and reflected-path interference 
of the uplink SI channel are both vanished if M and K satisfy 
M K and M is large enough. 

Proof: We consider the SI term in Eq. ® 


^W T G s Ax A = ^W T G s Ax A + ^W T G s Ax d . (17) 


For the ZF precoding and decoding case, by substituting Eq. 
© and Eq. © into Eq. (flTb . we obtain the first term in the 
RHS of Eq. (fT71 >. i.e., the direct-path SI 


^W t G s Ax d = a7lF ^{G H GY l G H G s G* (G^GT^d 

-1 


= VPd\ 


I M -1< 


G H G 

M 


k=l ^k 

/ s~iT \ -1 


G G S G* 

M 3 / 2 


G 1 G* 


x d . (18) 


All the factors will converge to a non-zero constant or constant 
matrices except the ( G R G J G ) when M is large enough. 
Rewrite G as the column form, then we have 


G H G S G* = [glg* 2 ,...,g* K YG s [gf9*2,-,9Ki (19) 


where g* is the 7th column of matrix G*, and 7 = 1,.... K. 
Note that G s satisfies the condition of Lemma |T| by assuming 
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all entries to be c max , namely, the worst case of direct-path 
SI. By using Lemma Q] we get 


g, G H g* 

M 3 / 2 


4 0, i,j = 1, K. 


( 20 ) 


Therefore, we have 


G H G S G* 

M 3 / 2 


0 KxK- 


( 21 ) 


Likewise, we can prove that the reflected-path SI term in Eq. 
(O is vanished. For the MRT/MRC processing case, we can 
multiply both sides of Eq. © by 1/M, then we obtain the 
result in a similar way. Thereby, with ZF or MRT/MRC linear 
processing, the SI will be vanished like the interference and 
noise in the TDD massive MIMO system. ■ 


B. The Uplink Decoded Signal 

With the large-scale antenna processing methods, the user 
data in the uplink received signal can be acquired when M is 
large enough. 

Proposition 1: In the shared-antenna full-duplex massive 
MU-MIMO system, if M and K satisfy M K and M is 
large enough, the uplink decoded signal for the /,'th user is 
given by 

rk -^4 s/Kxu,k, (22) 

for the ZF precoding and decoding, and 

y/p^xu , fcl (23) 

for the MRT precoding and MRC decoding. 

Proof: From Eq. © and Theorem |T] by using the 
asymptotic orthogonality property in (4), i.e., if alal^Pk 
and j-fgJg* "---> 0, i j, we obtain the results. ■ 

C. The Downlink Decoded Signal 

With the large-scale antenna processing methods, the down¬ 
link SI can also be suppressed. 

Proposition 2: In the shared-antenna full-duplex massive 
MU-MIMO system, if M and K satisfy M K and M is 
large enough, the downlink decoded signal for the kth user is 
given as 

—Vk -^4 y/pfx Aj k, (24) 

Pk 

where pk is given by Eq. CD. it is a linearly increasing 
function of \fM when M is large. 

Proof: The SI E^i (<=L + P'kqKq) and the noise term 
n A ,k in the RHS of Eq. (fill) are both Gaussian random 
variables with finite variance, and the number of UEs K is 
finite as well, if we multiply both side of Eq. CD by 4_, we 
get the result. ■ 

The result above looks fantastic. An intuitive explanation is, 
the received power at each user grows with the increase of the 
array scale M, while the downlink SI is bounded, therefore 
the user’s signal-to-SI ratio grows with the increase of M, 
namely, the SI can be suppressed by the large-scale antenna 
array. 



BS antenna array scale M 


Fig. 2. The full-duplex uplink (UL) and downlink (DL) SI or total interference 
plus noise versus the antenna array scale M. The normalized power of uplink 
and downlink are set as p u = 10 dB and p d = 13 dB, respectively. 


IV. Performance Results 

In the simulation, we neglect the insertion loss of the 
circulators and let the number of UEs K = 4. The large-scale 
fading /3k, (3, and /3' pq are set as 0.1, 0.8, and 0.7, respectively. 

Fig. ID simulates the uplink and downlink SI versus the 
antenna array scale M. The vertical axis represents the mean 
square (the power) of the SI and the total interference plus 
noise which are from Eq. © and Eq. CD - We can see that 
the SI and the total interference plus noise level get smaller 
with the increase of the array scale M, the SI level becomes 
lower with the decrease of the coefficient c r , or c’. As shown 
in Fig. [D the power of the SI converges to zero at a rate of 
1/M. It will decrease to the level 0 dB, i.e., the normalized 
white noise level, at M « 260 for the uplink with Cij = 0.5 
and M « 170 for the downlink with c' pq = 0.6. The total 
interference plus noise will decrease to 0 dB at M ss 430 
for the uplink (c l:] = 0.9) and M cs 220 for the downlink 
(c' pq = 0.7). Thereby, the SI can be suppressed effectively by 
the ZF and MRT/MRC methods when M is large enough. We 
also can see, the “large enough” M is sensitive to the direct- 
path SI coefficients c^j and c pq , thus the physical isolation 
or the duplex mode scheduling for the SI pre-suppression is 
necessary so that we can get a smaller but “large enough” M. 

V. Conclusion 

We proposed a shared-antenna full-duplex massive MU- 
MIMO model and discussed the uplink and the downlink self¬ 
interference suppression. This new model takes advantage of 
both the high spectral efficiency of the full-duplex system and 
the channel reciprocity of the TDD massive MIMO system. 
During the discussion of the SI suppression, we focus on the 
ZF and the MRT/MRC linear processing methods. A rigorous 
proof is given to demonstrate that the SI will be suppressed in 
the completely dependent uplink and downlink channels when 
the antenna array scale grows large enough. 
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